Entanglement versus chaos in disordered spin chains 
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We use a Heisenberg spin- 1/2 chain to investigate how chaos and localization may affect the 
entanglement of pairs of qubits. To measure how much entangled a pair is, we compute its concur- 
rence, which is then analyzed in the delocalized/localized and in the chaotic/non-chaotic regimes. 
Our results indicate that chaos reduces entanglement and that entanglement decreases in the region 
of strong localization. In the transition region from a chaotic to a non-chaotic regime localization 
increases entanglement. We also show that entanglement is larger for strongly interacting qubits 
(nearest neighbors) than for weakly interacting qubits (next and next-next neighbors). 
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I. INTRODUCTION 

Today it is well known that quantum entanglement, 
hereafter simply entanglement, is not just an issue of 
quantum mechanics that has yielded several discussions 
on the conceptual foundations of this theory 0, 0, . It 
is also a practical resource to develop new technologies. 
Entanglement enables us to implement some quantum 
algorithms which outperform their classical counterparts 
• It can also increase the amount of classical infor- 
mation transmitted in a quantum channel via the super 
dense coding protocol 0] and it allows the transmission 
of an unknown quantum state between two spatially sep- 
arated parties (quantum teleportation) 0, 01 ■ 

Since we can perform many useful tasks with entan- 
gled states, it is desirable to quantify the amount of en- 
tanglement these states have. There are at least three 
measures of entanglement which have reasonable phys- 
ical interpretations 0, 0, |n|. One of them, the En- 
tanglement of Formation (Ep), has a relatively simple 
analytical expression for bipartite mixed states [l^ . The 
Ep is a monotonically increasing function of a quantity 
called concurrence, which we adopt here as our measure 
of entanglement See section II for details. 

In this paper, we study how the entanglement of pairs 
of qubits is related to chaos and localization. To do so, we 
consider an one dimensional Heisenberg spin- 1/2 chain. 
The interest in such systems has recently increased con- 
siderably, for they model several proposed quantum com- 
puters [14( . Each spin in the chain corresponds to a qubit 
and the interaction between them is used to describe two- 
qubit gates. A single-particle excitation corresponds to 
a spin pointing up. In these systems, we call defect the 
site whose energy splitting differs from all the others. It 
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is obtained by applying a different magnetic field in the 
z direction to the chosen site. A disordered system is 
characterized by the presence of one or more defects. 

There have been recent attempts to relate localization 
with entanglement [l5| and also chaos with entanglement 
[l6l Il7j . Here we study at the same time the influence 
of localization and chaos on the entanglement. Some au- 
thors have shown that in a system of coupled quantum 
kicked tops G3 chaos increases entanglement. Here 
we show just the opposite, that chaos decreases entan- 
glement. Even though our system is very different from 
the one they considered, it is clear that a more careful 
analysis of the relation between chaos and entanglement 
is still needed. We claim that the entanglement of two 
qubits is not simply related to delocalization. It also is 
influenced by the chaoticity of the system. By chaos we 
mean how close the energy level spacing distribution of 
the system is to a Wigner-Dyson distribution. See next 
paragraph and Section III for details. 

The spin chain considered in this paper is a per- 
fect system to study the relations between entanglement 
and localization and between entanglement and chaos, 
because it has different regions of interest. Concur- 
rences of pairs of qubits obtained in different regimes 
can then be compared. The transition from integra- 
bility to non-integrability in these systems depends on 
the defects [2(J. It is known that the energy level 
spacing distribution of an integrable system is Poisso- 
nian, Pp(s) = exp(— s), while the level statistics of a 
chaotic system is given by the Wigner-Dyson distribu- 
tion, Pwd{s) — (tts/2) exp(— 7rs 2 /4) |2j|. In the absence 
of defects, that is, when we have an ideal chain, the sys- 
tem is integrable and it can be solved with the Bethe 
ansatz [22^ . Its level distribution is therefore Poissonian. 
As random on-site magnetic fields are turned on and their 
mean-square amplitude starts increasing the system un- 
dergoes a transition and becomes chaotic. In this sce- 
nario we obtain a Wigner-Dyson distribution for the level 
spacings. By further increasing the mean-square ampli- 
tude, localization eventually takes place and the distri- 
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bution becomes Poissonian again. 

When the states of the system become localized, the 
excitations in the chain are restricted to finite regions 
of space. Defects even more different in energy finally 
leads to strong localization, by which we mean that 
each eigenvector of the system becomes very close to a 
non-interacting multi-excitation state \<j>j) — \a±, ...,ocl), 
where otk = 0, 1 indicates a spin down or up, respectively, 
and L is the total number of sites. Each excitation is 
then restricted to practically just one site of the chain. 
In quantum computation, these are the states we want to 
operate with and they are called quantum registers |23| . 
On the other hand, the eigenfunctions of the system in 
the absence of defects and especially the eigenfunctions in 
the chaotic regime are very much spread over the chain 
and they correspond to linear superpositions of several 
quantum registers. The system in this case is said to be 
delocalized. 

The Hamiltonian that describes our system is the fol- 
lowing, 

n— 1 n— 1 

L-l 

+ + (i) 

n— 1 

We consider only nearest neighbor interaction, h = 1, 
and a x , a v , and a z are Pauli matrices. Each site n is 
subjected to a magnetic field in the z direction, giving 
the energy splitting h n . There are L sites and we deal 
with an open chain, which is, in our opinion, a more 
realistic model for quantum computers. The system is 
isotropic (anisotropic) when the coupling constant Jz for 
the diagonal Ising interaction <7*<7^ +1 is equal (different) 
to the coupling constant Jxy for the XY-type interac- 
tion <T'n<rf l+1 + ^n a n+i- This last term is responsible for 
delocalizing the system, because it propagates the exci- 
tations through the chain. It is usually referred to as the 
hopping term. 

To analyze how localization affects the entanglement 
of a pair of qubits, we compare the number of principal 
components for the whole system with the concurrence 
for a chosen pair. The definitions for these two quanti- 
ties are presented in the next section. There we restrict 
ourselves to two qubits and show a brief analytical dis- 
cussion. As expected, entanglement tends to disappear 
in a strongly localized system. On the other hand, it be- 
comes maximal if the level spacing of the two qubits is 
the same and it is also very large when the interaction 
between them is much larger than the difference between 
the two level spacings. 

In section III, we show our numerical results for a chain 
with several qubits and several excitations. First, we 
discuss the case where two qubits have the same level 
spacings, which are in turn very different from the other 
qubits. If the Ising interaction does not exist and only the 
hopping term is present, we have a situation very similar 



to the case of only one excitation discussed in section II 
and maximal entanglement can be obtained. However, 
if the Ising term is present, the entanglement for nearest 
neighbor qubits becomes larger than for the other pairs, 
indicating that many-body effects can strongly affect en- 
tanglement between distant neighbors. 

In section III, we also consider qubits with random 
level spacings, which allows us to study the behavior of 
entanglement with chaos and the behavior of entangle- 
ment with localization when several interacting excita- 
tions are present. In the integrable but delocalized re- 
gion of an ideal chain, entanglement is large. It then de- 
creases when the system becomes chaotic and even more 
delocalized. However, in the transition region between 
chaos and strong localization, where the system is actu- 
ally localizing, entanglement increases again. It is only 
when the system gets close to strong localization that 
entanglement decreases. The relation between localiza- 
tion and entanglement is not so simple as one might have 
expected. But it is clear from our analysis that chaos di- 
minishes entanglement. We should add, however, that 
for Hamiltonian Q nearest neighbor interacting qubits 
are able to keep their concurrence reasonably large even 
in a chaotic region. 

Concluding remarks are presented in section IV. 

II. ANALYTICAL DISCUSSION: TWO QUBITS 

In order to study quantitatively the entanglement be- 
tween two pairs of qubits we should calculate their en- 
tanglement of formation (Ep). Given the density matrix 
P12 that describes our pair of qubits, Ep is the average 
entanglement of the pure states of the decomposition of 
P12, minimized over all possible decompositions: 

E F (p 12 ) = min^p^(^), (2) 

i 

where J2iPi = 1> < Pi < 1) and p 12 = J^iPi IV'i) (V'il- 
Here E(ijj) is the von Neumann entropy of either of the 
two qubits |24|. Wootters et al 01 have shown that, for a 
pair of qubits, Ep is a monotonically increasing function 
of the concurrence, which one can prove to be an entan- 
glement monotone. Since the concurrence is mathemat- 
ically simpler to deal with than Ep, we concentrate our 
efforts calculating the concurrence to study the amount 
of entanglement between two qubits. The concurrence 
between them is [l2|: 

C12 = max{Ai - A 2 - A 3 - A 4 , 0}, (3) 

where Ai, A2, A3, and A4 are the square roots of the eigen- 
values, in decreasing order, of the matrix R = pnpn- 
Here hi is the time reversed matrix given by 

PX1 = ^ CT y ) Pl 2 (<7y (g) <Jy) . (4) 

The symbol p* means complex conjugation of the matrix 
p in the basis {|11) , 1 10) , |01) , |00)}. 
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The above procedure to calculate the concurrence of 
two qubits is used in the next section, where we deal with 
several qubits. There, we trace over the qubits we are not 
interested in and study the reduced density matrix of the 
two chosen ones. For a pure bipartite system, which is 
the case for the eigenvectors of a two spin chain, the 
concurrence is simply given by: 



C12 =2\ad- bc\ , 



(5) 



where |V>) 12 = a \U)+b |10)+c \01)+d |00) is the bipartite 
pure state. We see that when the concurrence is zero the 
two qubits are not entangled and when it is 1 they are 
maximally entangled (EPR-Bell states). 

For an isotropic spin chain with two sites the Hamil- 
tonian can be written as: 



H 



\ (h x a\ + h 2 a z 2 ) + J - Ka| + + of of ) . (6) 



One can easily verify that [H, erf +<rf] = 0, i. e., the total 
angular momentum in the ^-direction is conserved. This 
implies that states with different number of excitations 
do not couple. This is evident by looking at the matrix 
form of H: 



H = 



V 







(7) 



where £ = hi + h<x and A = h\—h%. The above Hamilto- 
nian is block diagonal, which means that |11) and 1 00) are 
non-entangled eigenvectors of H . As we are interested in 
entanglement, we restrict our analysis to the other two 
eigenvectors. These two eigenvectors are obtained via the 
diagonalization of the above 2x2 block matrix. After a 
straightforward calculation we obtain the two remaining 
eigenvectors and eigenvalues: 



\E±) = 



1 



2(J 2 + A 2 )t2AVJ 2 + A 2 

j|io>- (a t Vj 2 + aA |01)1 , (8) 



For any other value of A we cannot achieve maximal 
entangled states, showing that the appearance of a defect 
(hi 7^ /12) in this two-site chain reduces the amount of 
entanglement. But, even if A ^ 0, when J » A wc 
get large concurrences. This means that in the regime of 
strong coupling we can still have large entanglement. 

To investigate the relationship between entanglement 
and localization we compute the number of principal 
components, N pc . This is a quantity often used to de- 
termine how much spread each eigenstate of the system 
is. The eigenvectors of the Hamiltonian (JIJ are written 
as linear superpositions of the quantum registers \<fii), 
i.e. \ipj) = J2i=i n a i \4>i)i where N is the total number 
of eigenstates. The N pc for the eigenvector j is defined 
as H3 



N 3 

pc 



1/ 



N 



(11) 



A delocalized system has large N pc , while a strongly lo- 
calized system has N pc very close to 1. 

Applying Eq. ltTT|) to the eigenvectors \E±) we get the 
following expression for N pc : 



JV* = - 



Ci/2 



(12) 



Eq. (|12|) shows that N pc for these eigenvectors is an in- 
creasing function of the concurrence, which implies that, 
for this particular two-qubit system, the more it is delo- 
calized the more it is entangled. When we have C± — 1 
(maximum entanglement), N^ c — 2. When C± = 
(no entanglement) we get A/^ = 1 (strong localization). 
However, this association between derealization and en- 
tanglement is not always valid and this will become clear 
in the next section, where we deal with several qubits. 
For the moment, it is worthy noting that for a gen- 
eral superposition of the four quantum registers we may 
have large N pc but no entanglement at all. Even though 
the state = y/T/l (|11) + 1 10) + |01) + |00>) is not an 
eigenstate of our Hamiltonian, it can still be used as a 
good example of this situation, for it has N pc = 4 but 
C v , = 0. 



± 4 2 



(9) 



Using Eq. JSJ) we get the following expression for the 
concurrence of the eigenvectors \E±): 



C± 



v/l + A 2 /J 2 



(10) 



Eq. i|rU|) clearly shows that C± is a decreasing function 
of A and an increasing function of the coupling constant 
J. When A = 0, that is when the two qubits have the 
same level spacing, we obtain maximal entangled states 
whether we are in the weak or strong coupling regime. 



III. NUMERICAL RESULTS: SEVERAL 
EXCITATIONS 

In this section we study how chaos and localization 
may affect the entanglement of a pair of qubits when 
several qubits and several excitations are considered. De- 
pending on the defects, our disordered system may be- 
come chaotic. To determine if a system is chaotic or not 
it is usual to compute its energy level spacing distribu- 
tion. Notice that in the model described by Hamiltonian 
JIJ, the z component of the total spin J2t=i ^ s con " 
served, so states with different number of excitations are 
not coupled. We therefore analyze sectors with the same 
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number of excitations. Since we are interested in com- 
paring chaos and entanglement, we focus on the sector 
with the largest number of states, that is, the sector with 
L/2 excitations, because this is the region where chaos 
should set in first |23| . 

In a very large system, the boundary conditions have 
no effects, but numerical calculations are limited to a fi- 
nite number of sites. A chain with free boundaries (open 
chain) and defects only on the edges is known to be in- 
tegrable [2(|. For our numerical calculations we choose 
an open chain with defects of values — J/2 on the edges. 
Such values should diminish border effects. The Hamil- 
tonian becomes now 



H = H 



Jz (1 



•Mi 



(13) 



We work with L = 12 sites and 6 excitations, which gives 
us 12!/(6!6!) = 924 states. 

Before analyzing how chaos can affect entanglement 
and motivated by the previous section, we first check if 
two qubits with the same level spacing, but now in a 
chain with several excitations, can also lead to maximal 
entanglement. Suppose that two selected qubits have 
level spacing h + d, while all the others have level spacing 
h. If we turn off the Ising interaction (i.e., if Jz = 
in Eqs. JQ) and l|13l) h and if d ^ Jxy, the two selected 
qubits indeed lead to maximal entanglement (except if 
one qubit is on the edge of the chain, in which case border 
effects are noticed). This can be seen from the top left of 
Fig. QJ, where circles indicate that the chosen pairs are 
the nearest neighbor qubits (n, n + 1), squares give next- 
nearest neighbors (n, rt + 2), and triangles give next-next- 
nearest neighbors (n, n + 3). In the figure we set Jxy = 1 
and d = 100. The concurrence is obtained by tracing 
over the qubits we are not interested in. We show the 
maximum concurrence Cmax obtained among the 924 
eigenstates for each selected pair. When Jz — 0, even 
though we have many excitations, the situation resembles 
the case of one single excitation, where only the hopping 
term is present. Similar to this discussion on 

how to tune chosen qubits in order to find maximally 
entangled states was done in p7|. 

However, when the Ising interaction is on, the results 
can change considerably. Many-body effects now play 
an important role. For the case of an isotropic chain 
(Jz = Jxy = J) shown on the top right of Fig. QJ, only 
nearest neighbors maintain large concurrence, though not 
maximal. A particular case of anisotropy (Jz > Jxy), 
on the other hand, tends to increase the concurrence 
even for next-nearest neighbors, as can be noticed by 
comparing the graphs in the middle of Fig. QJ. On the 
left Jz = IOJxy and on the right Jz = IOOJxy- But 
the anisotropy has unexpected effects on the next-next- 
nearest neighbors. In the especial case of Jz = d = 
100 Jxy shown in the middle right of Fig. JIJ, their con- 
currences are also increased, but this is not the case for 
large anisotropy where Jz =/= d. In this situation, the 
concurrence of some pairs can become very close to 1, 



while others may decrease a lot, as can be seen in the 
graphs at the bottom of Fig. QJ, where arbitrary large 
values of Jz were chosen (Jz = d = 159Jxy on the left 
and Jz = d = 327 Jxy on the right). Since the main goal 
of this paper, however, is to analyze the relation between 
chaos and entanglement, we will let for a future publica- 
tion a more careful study of the many-body effects due 
to the relative strength of the Ising (Jz) and hopping 
(Jxy) terms in the Hamiltonian JIJ. 
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FIG. 1: (Color online) We show the maximum concurrence for 
several pairs of qubits with the same level spacing. Circles in- 
dicate that the selected pairs correspond to the nearest neigh- 
bor qubits (n, n + 1), squares give the next-nearest neighbors 
(n, n + 2) , and triangles correspond to the next-next-nearest 
neighbors (n,n + 3). We set d — 100 and Jxy = 1. Top 
left: Jz = 0, top right: Jz = Jxy, middle left: Jz = IOJxy, 
middle right: Jz = 100 Jxy, bottom left: Jz ~ 159 Jxy, and 
bottom right: Jz = 327 Jxy- 

We now proceed to study the relation between local- 
ization and entanglement and also between chaos and 
entanglement in our many-body system. We consider an 
isotropic chain (Jz = Jxy = J) and set J = 1. The 
level spacings of the qubits are obtained with random 
magnetic fields along the z direction. They are given 
by h n = h + d n , where cf n 's are uncorrelated random 
numbers with a Gaussian distribution: (d n ) — and 
(dndm) — d 5 n ^ m . 

According to the integrability and localization of the 
system, we can identify different regions. This is shown 
in Fig. J2J. We calculate the average of the number 
of principal components N pc for the 924 eigenstates as 
a function of d. This is done using a sequence of 12 
Gaussian random numbers, which give the random level 
spacings for the qubits. This procedure is then repeated 
for 20 different sequences. At the top of Fig. J2J we 
show (N pc ), where ( ) corresponds to the average over 
these 20 different sequences. We decided to work with 
20 sequences, because this is enough to give an idea of 
the general behavior of the system. We compared the 
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results for the concurrences of some pairs for more se- 
quences and the behavior was still very similar, there- 
fore justifying the use of only 20 sequences. At the bot- 
tom of Fig. J2J we show a quantity used to measure how 
much chaotic the system is. This parameter is defined 
as V = f°°[P(s) - P WD (s)]ds/ f°°[Pp(s) ~ P WD (s)]ds, 
where so = 0.4729... is the intersection point of Pp(s) 
and Pwd (23l l2l| and P(s) is the energy level spacing 
distribution for the system under study. A regular sys- 
tem has 77 = 1 and a chaotic system has r\ — 0. Here 
again, ( ) indicates the average over 20 sequences of ran- 
dom numbers. 

When d = the system is integrable, but delocalizcd. 
We have Poisson distribution (77 ~ 1), but large (N pc \. 
As d increases the system becomes chaotic and even more 
delocalized. For < d < 0.2, we move toward a Wigner- 
Dyson distribution (77 tends to 0) and (A pc ) becomes 
even larger. However, as we further increase d the level 
spacing distribution approaches again a Poissonian distri- 
bution and (N pc ) decreases. This transition region cor- 
responds to 0.2 < <i < 2. As d becomes much larger than 
the interaction strength J, the system becomes strongly 
localized, the distribution is again clearly Poissonian and 
(N pc ) gets very close to 1. 
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d/J 

FIG. 2: Top: Dependence of (N pc } on d/J. The over line 
corresponds to an average over the 924 eigenstates for each 
d/J and ( ) comes from the next average over 20 different 
sequences of random numbers. Bottom: Dependence of {77) 
on d/J. Here again, { } indicates an average over 20 different 
sequences of f 2 Gaussian random numbers. We set J = 1. 

We compare (N pc ) and (77) with the concurrence for 
pairs of nearest neighbors (Fig. ©), pairs of next-nearest 
neighbors (top of Fig. 10}), and pairs of next-next-nearest 
neighbors (bottom of Fig. (0}). For each pair we compute 
the maximum concurrence C m ax among the 924 eigen- 
states. This is averaged over 20 different sequences of 
random numbers. We also compute the average concur- 



rence C for each chosen pair, which is again averaged over 
20 different sequences of random numbers. The average 
concurrences for all pairs have a very similar behavior. 
Because of this and also because the values for the ave- 
rage concurrences become very small for distant pairs, 
we just show here the average concurrences for nearest 
neighbor qubits (bottom of Fig. @). 
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d/J 

FIG. 3: (Color online) The top gives the dependence of 
(Cmax) on d/J and the bottom gives (C) for nearest neighbor 
qubits (N). In both: solid and black line represents the pair of 
qubits 1-2, dotted and red line gives the pair 3-4, dashed and 
blue gives the pair 6-7, and dot-dashed and green represents 
the qubits 10- If . 

By using Figs. J2J), 10 , an d (QJ we can analyze what 
happens to entanglement in the distinct regimes of the 
system. 

First, let us look at the bottom of Fig. ©. By com- 
paring the average concurrences (C) for d — and 
< d < 0.2, the region where chaos and derealization 
increase, we see that the concurrences slightly decrease. 
This suggests that in such systems, with nearest neigh- 
bor interactions, chaos may contribute to a decrease of 
the average entanglement of the qubits. In the transi- 
tion zone, 0.2 < d < 2, which is the region where the 
system becomes less chaotic and more localized, we see 
a rapid increase of the concurrences. This is very inter- 
esting, because despite of being in the process of local- 
ization the average concurrences of the qubits increase. 
Finally for d > 2, the region where the system is strongly 
localized and clearly regular (non-chaotic), the average 
concurrences decrease, as it should. The effect of local- 
ization on the average entanglement, therefore, depends 
on how far we are from the chaotic region. As explained 
above, localization increases entanglement when the sys- 
tem is moving from non-integrability (chaos) to integra- 
bility, but it decreases entanglement when the system is 
already strongly localized. 

The top of Fig. J3J shows that in the region of chaos 
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FIG. 4: (Color online) Top: Maximum concurrence for next- 
nearest neighbor qubits (NN). Solid and black line gives the 
pair 1-3, dotted and red is for the qubits 2-4, dashed and 
blue represents 3-5, and dot-dashed and green gives the pair 
4-6. Bottom: Maximum concurrence for next-next-nearest 
neighbor qubits (NNN). Solid and black line gives the pair 
2-5, dotted and red is for the qubits 4-7, dashed and blue 
represents 5-8, and dot-dashed and green gives the pair 7-10. 

(d ~ 0.2) the maximum concurrences for nearest neigh- 
bor qubits are larger than (C max ) for next and much 
larger for next-next nearest neighbors (Fig. J3J). This is 
explained noticing that for directly coupled pairs (near- 
est neighbors), and therefore pairs with strong interac- 
tion, the effects of chaos are not so drastic and they keep 
reasonably large maximum concurrences. But again we 
verify that chaos diminishes the entanglement of any pair 
of qubits. We also see that in the region where localiza- 
tion is becoming strong (d > 2) we still have reasonable 
values for (C max ) for the nearest neighbor qubits. This 
contrasts with the behavior of (C max ) for the next and 
next-next nearest neighbor qubits, where we see a faster 
falling (C max ). 

From the observations of the previous paragraph we 
conclude that the interaction between two qubits coun- 
terbalances the destruction of their entanglement caused 
by chaos and also by strong localization. This inter- 
pretation is reinforced if we note that the Hamiltonian 
considered here has only nearest neighbor interaction 



terms and, as it was numerically shown, nearest neigh- 
bors maintain larger entanglement if compared to next 
and next-next neighbors, which are shown to be more 
susceptible to the effects of chaos and localization. 



IV. CONCLUSION 

We studied how entanglement is related to chaos and 
localization in one dimensional Heisenberg spin chains. 
For this purpose we used an isotropic chain with external 
random magnetic fields (chain with defects) to obtain a 
chaotic regime. 

We showed that chaos is responsible for a decrease 
in the entanglement (concurrence) of the nearest, next 
and next-next-nearest neighbor qubits. However, near- 
est neighbor qubits are less sensitive to chaos than the 
next and next-next-nearest neighbors and keep relative 
large entanglement in the chaotic region. 

The relation between entanglement and localization is 
rather subtle. We found that in a strongly localized sys- 
tem entanglement decreases. Nonetheless, in the transi- 
tion region from non-integrability (chaos) to integrability 
we observed that an increase in localization causes an in- 
crease in entanglement. 

Finally, it was also shown that when two qubits have 
the same external magnetic field and this field is very 
large compared with the fields applied to the other sites, 
they lead to maximally entangled states if the Ising inter- 
action is not present. However, this behavior is affected 
by the relative strength of the Ising (Jz) and hopping 
{Jxy) terms in the Hamiltonian. Large anisotropy leads 
to large entanglement for nearest and next-nearest neigh- 
bor qubits, but the behavior of the next-next-nearest 
neighbors does not follow this or any simple trend. 
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